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Abstract. We revisit the problem of computing the spreading and covering numbers. 
These numbers were first introduced by Geramita, Gregory, and Roberts to study the 
Ideal Generation Conjecture, a conjecture concerning sets of points in a projective space. 
Both the spreading number and covering number can be reformulated as a graph theory 
question. We describe some new approaches to computing or bounding these numbers. 
In addition, we show a surprising connection to some of the spreading numbers to the 
sequence A053307 in the On-Line Encyclopedia of Integer Sequences. 



1. Introduction 

Let R = k[xi, . . . , Xn] be a polynomial ring over a field k. For any non- negative integer 
d, let Mrf be the set of all monomials of degree d in R. For any subset W C M^, let 

RiW = {xiTn \ m & W and I < i < n}, 

that is, RiW is the subset of M^+i formed by multiplying each monomial of W by each of 
the variables of R. For any W C Md, we always have |-RiVr| < n\W\ and RiW C Md+i- 

We are interested in finding subsets W where either |-RiVr| = n\W\ or RiW = M^+i. 
In particular, we define the spreading number to be 

an{d) = max{ \W\ : W C Md and \RiW\ = n\W\}. 

The terminology is derived from the fact that the elements of RiW are "spread" out in 
Md+i- Similarly, the covering number is defined to be 

p^{d+l) = mm{\W\ : W C Md and RiW = Md+i} . 

In this case the elements of RiW "cover" the elements of Md+i- The goal of this paper is 
to introduce some new approaches to computing and bounding these numbers. 

The spreading and covering numbers were introduced by Geramita, Gregory, and 
Roberts to study the Ideal Generation Conjecture for a set of generic points in P" (see 
[3l Theorem 4.7]). Geramita, et al. gave exact values for an{d) for all d when n = 3 
or 4, and some scattered results and bounds on other values. Curtis [2] later found a 
formula for psid) for all d and an improved lower bound on Pi{d). Using techniques from 
linear programming, Hulett and Will [4J improved these lower bounds on Pi{d). Carlini, 
Ha, and the second author pj later reformulated the problem to show it is equivalent to 
computing the dimension of some abstract simplicial complex. 
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The computation of these invariants has proven to be quite challenging. Computing 
an{d) or pn{d + 1) by brute force, even for small n, quickly proves infeasible due to the 
number of subsets one must check. For example, to compute 05(5), we have \Md\ = 
_ j^2g^ so there are 2^^^ ^ 8.5 x 10^"^ subsets that need to checked using a 
naive method. We have therefore focused on introducing algorithms that bound an{d) 
and Pn{d+ 1). Unfortunately, as n and d increase, so does the distance between our two 
bounds. Given the difficulty of computing these invariants, we have therefore set ourselves 
a goal of bounding as many of the spreading and cover numbers in a neighbourhood of a 
size of no more than 100. 

To find these bounds, we make use of a construction of Geramita, et. al; in particular, 
we make a graph Sn{d) where the vertices of this graph correspond to the elements of 
Md- The values of and pn{d + 1) are then related to graph theoretic invariants of 

Sn{d). In Section we make this connection explicit, and reframe our problem in terms 
of Sn{d). We also highlight the symmetry of Sn{d), which we exploit in our algorithms. 

In Section [31 we focus on an algebraic method to compute an upper bound on Q;„((i). 
We associate to Sn{d) its edge ideal I{Sn{d)). Computing an{d) then reduces to com- 
puting the KruU dimension of T / I{Sn{d)) where T is some suitable polynomial ring. We 
construct suitable linear forms Li, . . . , L^, to make the ring T / [I[Sn{d)) + (Li, . . . , Ls)). 
The dimension of this new ring allows us to give an upper bound 

In Section m we describe two greedy algorithms. The first algorithm gives a lower bound 
for while the second gives an upper bound for pn{d+l). Both exploit the symmetry 

of the graph Sn{d). We also include refinements to both algorithms so that they can be 
adapted to compute and Pn{d + 1) exactly. We summarize our output in Section [5] 

where we compare output to the bounds found in previous papers. 

In Section E] we prove that the sequence 04 (rf), one of the few cases for which we have 
a formula, equals the sequence A053307 found on the On-Line Encyclopedia of Integer 
Sequences fl\. This sequence counts the number of 2 x 2 matrices with non- negative 
entries that sum to d up to permutation of rows and columns. One is lead to ask whether 
other spreading or covering numbers are related to interesting integer sequences. 

As a final note, all of the code that we used will be available on the authors' website^. 

Acknowledgements. This work was made possible by the facilities of the Shared Hi- 
erarchical Academic Research Computing Network (SHARCNET:www.sharcnet.ca) and 
Compute/Calcul Canada. We used Macaulay 2 [6] to carry out our calculations. The first 
author was supported by an NSERC USRA during the summers of 2010 and 2011. The 
second author was supported by an NSERC Discovery Grant. 

2. Preliminaries 

We recall the relevant definitions and preliminary results that we will need for the 
remainder of the paper. As in the introduction, we let R = k[xi, . . . ,Xn] denote the 
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polynomial ring over a field k, and Md will denote the monomials of degree d in R. We 
first translate our problem of computing and pn{d+l) into a graph theory problem. 

Construction 2.1. Fix positive integers n and d. Let Sn{d) denote the graph whose 
vertex set is the set of monomials in R = k[xi, . . . ,x„], and two vertices mi,mj are 
adjacent if and only if deg ( 1cm (mj, mj)) = d + 1. 

Note that we will abuse notation and use to mean both the vertices of Sn{d) and 
the set of monomials of degree d in R = k[xi, . . . , Xn]- 

Example 2.1. If we consider 5*3(3), then xlx2 and xfxs are adjacent but xfa;2 and x^x^ 
are not. The graph S^lS) is given in Figure 1. 




O , >-» 2 ,-->-» ^ 0''2 /y>3 



Figure 1. ^3(3) 

The spreading and covering numbers are related to the maximum independent set and 
a special type of minimum clique covering, respectively, of Sn{d). 

Definition 2.2. A subset V C is an independent set if any two distinct elements 
of V are not adjacent; is a maximal independent set if it is not properly contained 
in any larger independent set. 

Definition 2.3. A subset of in which any two vertices are adjacent is called a clique. 
A clique is maximal if it is not properly contained in any larger clique. If Ci,...,Ct 
are cliques, we say they form a clique cover of Sn{d) if Ci U ■ ■ ■ U = Md. For any 
monomial m of degree d — 1, an upward clique is the clique consisting of the vertices 
mxi e Md for all Xi G {xi, . . . , x^}. 

As shown in [3], an{d) and Pn{d + 1) are equivalent to an invariant of Sn{d): 

Lemma 2.4. With the notation as above 

(z) an{d) is the cardinality of the largest maximal independent set of Sn{d). 
{a) pn{d + 1) is the minimum cardinality of an upward clique cover of the vertices of 

Sn{d+l). 
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Example 2.5. The graph 5*3(3) in Example 12. II has 03(8) = 4 because {xf , x^, x^, X1X2X3} 
forms a maximal independent set. For the same graph, we have ^3(2 + 1) = 4 be- 
cause Ci = {xl,xlx2,xlx3}, C2 = {xl,Xixl,X3xl}, C3 = Xi^l, Xaxl}, and C4 = 
{x'lx2, X1X2, 3^13^23^3} form a minimal upward clique cover of 5*3(2 + 1). 

Finding the size of the largest maximal independent set or a minimum clique cover of 
a graph are both NP-hard problems. This explains, in part, why computing an{d) and 
pn{d + 1) is so difficult. However, the symmetry of 5*„((i) offers hope for using specialized 
techniques to improve our bounds on a;„((i) and Pn{d + 1). 

Let Sym(n) denote the symmetric group on the set {1, 2, . . . , ra}. For any x** = 
2^1^ ■ ■ ■Xn"'" e and a G Sym(n), let (t(x'*) denote the monomial obtained by permut- 
ing the indices 1, . . . ,n according to the permutation a. This operation preserves many 
properties of sets of vertices; e.g., independent sets and clique covers are both unaffected. 

Lemma 2.6. IfW = {mj^, . . . is an independent set on Sn{d), then the set 

a{W) = {cT{mj) I TTij e W} 
is also an independent set of Sn{d) for any a G Sym(n). 

Proof Since W is an independent set, for any two mj,mk G W, deg ( lcm(mj, m^)) 7^ d+1. 
Applying a to each element in W results only in an exchange of indices; the exponents of 
each indeterminate, and thus the degree of the monomial, remain unchanged. Therefore, 
deg ( lcm(cr(mj), cr(mfc))) ^ d+1, so a{mj) and (rijUk) are not adjacent. □ 

The proof that clique covers are preserved is similar. We will elaborate further on how 
the symmetry of Sn{d) comes into play in Sections [3] and |H For now we present a recursive 
lower bound for an{d) that uses another fact about the structure of Sn{d). 

Theorem 2.7. For all n > 2,d > 3, an{d) > an{d - 2) + an-i{d) . 

Proof. Consider all the vertices of Sn{d) divisible by xf, i.e., those vertices which are 
labelled with a monomial of the form x^^ ■ ■ ■ with ai > 2. The induced graph on these 
vertices is isomorphic to Sn{d — 2). Also, consider all the vertices which are not divisible 
by Xi. The induced graph on these vertices is isomorphic to Sn-i{d). 

The graph Sn{d) consists of these subgraphs, plus all the vertices whose corresponding 
monomials have the form xiXg^ ■ ■ ■ x^", i.e., the degree of xi is one. These extra monomials 
form a "buffer" between the two subgraphs, that is, no vertex in the subgraph 5*„((i — 2) 
is adjacent to a vertex in the subgraph Sn-i{d). Thus, the union of an independent set in 
Sn{d — 2) and an independent set in Sn-i{d) will be an independent set in Sn{d). From 
this observation we get an{d) > dnid — 2) + an-i{d). □ 

We will compare various lower bounds for an{d) in Section [51 In general, the recursive 
lower bound is not as tight as the lower bound found in [3]. 
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3. Upper Bound for the Spreading Number 

We will describe how to find an upper bound on by using some techniques from 

commutative algebra. In particular, we first show that can be encoded as the 

dimension of particular ring and then describe how to bound this dimension. 

We first take our graph theory problem of Lemma 12.41 and translate it into an new 
question. For any finite simple graph G with vertex set Vq = {zi, . . . , Zn} and edge set 
Egi we can associate to G a quadratic square-free monomial ideal called the edge ideal. 
In particular, given G, the edge ideal of G is then 

I{G) = {ziZj I {zi, Zj} G Eg) C T = k[zi, ...,Zn]. 

Understanding how the graph theory invariants of G are encoded into the invariants of 
T/I{G) is an ongoing area of research (e.g., see p]). In particular, it is known that 

dimT//(G) = a{G), 

where a{G) is the independence number of G, that is, is the cardinality of the largest 
independent set. When G = Sn{d), it follows by Lemma that a{Sn{d)) = an{d), the 
spreading number. We have thus proved: 

Lemma 3.1. Let I{Sn{d)) = {zj^Zmj \ {tTii,mj} is an edge of Sn{d)) in the ring T = 
k[zm I m G Md]. Then 

Cin{d) = dimT/ I{Sn{d)). 

We therefore want to compute dimT/I{Sn{d)). Computer algebra systems, such as 
CoCoA or Macaulay2, normally tackle this problem by computing the associated Hilbert- 
Poincare series. This approach improves upon naive methods, but even for small cases, 
many computer algebra systems cannot compute this dimension. 

Remark 3.2. Computing an{d) by computing dimT/ I{Sn{d)) was first described in 
[Ij. Our initial hope was that improved computer hardware in the intervening ten years 
would allow us compute new values of an{d). This hope ended up being too optimistic. 
For example, computing 05(5) in Macaulay 2 requires us to compute the Hilbert-Poincare 
series of an ideal in a ring with (^+5^^) = 126 variables. 

Instead of using the Hilbert-Poincare series approach to computing the dimension of a 
ring, we wish to make use of a system of parameters. 

Definition 3.3. Let / C i? = k[xi, . . . homogeneous ideal. A partial homoge- 

neous system of parameters (hsop) of R/I is a finite sequence of homogeneous forms 
Fi, . . . ,Ft with deg Fi > for each i such that 

dim(i?/(/, Fi, . . . , Ft)) = dim R/I - t. 

A partial homogeneous system of parameters is a (complete) homogeneous system 
parameters if t = dim R/L 

If we could find a hsop for T/ I{Sn{d)), then we would be able to compute dimT/I{Sn{d)). 
However, the standard algorithms for computing a hsop (see [5]) require knowing dim R/I 
in order to compute its hsop. We can still bound an{d) using the following lemma. 
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Lemma 3.4. Let Li, . . . , Lt be any t linear forms ofT. Then 

dimT/{I{Snid)), L^,...,Lt)+t> dimT/I{Sn{d)) = 
Moreover, we have equality if Li, . . . ,Lt forms a partial hsop. 

Proof. The second part of the proof follows directly from the definition of a partial hsop. 
The first part of the proof follows from the more general fact that for any homogeneous 
ideal / in T and linear form L G T, then dimT/(/, L) > dim T/J — 1. □ 

This lemma forms the basis of our algorithm to compute an upper bound on Q;„((i). 
We want to judiciously pick a set of linear forms C = {Li, . . . ,Lt} so that computing 
dimT/ I{Sn{d), Li, . . . ,Lt) is easier than computing dimT//(5'„((i)). At the same time, 
we want to pick our Lj's so that a large subset of C forms a partial hsop. 

For the moment, let us suppose that we know how to pick £; we can then present the 
pseudo-code for this approach: 

Algorithm 3.5. Compute an upper bound for a;„((i) 

Input: n,d — The number of variables and degree of monomials, respectively. 
Output: An upper bound for an{d) 

Step 1 Let T = ^[2:^1, . . . , ^^mj where = {mi,...,m^} is the set of {^~^'^~^) = ^ 

monomials of degree d, and let I{Sn{d)) be the edge ideal of Sn{d) in T. 
Step 2 Pick a suitable choice of linear forms C = {Li, . . . , Lt}. 
Step 3 Compute dimT/J where J = L{Sn{d)) + (L | L G C). 
Step 4 Return dimT/J + 

We still need to explain how to pick C in Step 2 of Algorithm 13.51 We describe three 
strategies and comment on their strengths and weakness. 

3.1. Random Linear Forms. Fix integers n and d, and assume that some oracle has 
given a crude lower bound on (e.g., [3] showed that Vn{d)/n < an{d) where Vn{d) = 

{"^~^d~^) ~ Pick t to be equal to this lower bound. Form a t x Vn{d) matrix C 

with random entries taken from the field k, and let C be the following t linear equations: 



C 



We then make use of the following theorem of Kemper: 

Theorem 3.6 ([5l Proposition 1]). Suppose that A = R/L has dimension n, and let 
di, . . . ,dt G N>o. Then the set 

{{F^,...,F^)■. dim(i?/(/, F^, . . . , F^)) =n-t} 

forms a Zariski open subset of x ■ ■ ■ x A^t ■ 

Our linear forms C, for a suitable choice of matrix C, belong to the Zariski open 
subset Ai X ■ ■ ■ X A^. In other words, we expect most, if not all, of the elements of 
L to form a partial hsop. Indeed, our computer experiments on small values of n and 
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d appear to support this statement. However, the computation of dimT/J with J = 
I{Sn{d)) + (L I L G £) in Step 3 of our algorithm appears just as difficult as computing 
dimT/ I{Sn{d)). We had to abandon this approach because it was not practical. 

3.2. Using the symmetry. The second construction of C makes use of the symmetry 
of the graph Sn{d). Consider the cycle a = {12- ■■n — 1) G Sym(r;,). For each vertex 
m G Md of the graph Sn{d), define the linear form: 

L{m) = Zm + Z„{rn) H h Z„r.-i^^) G T. 

Here, z^, is the variable of T = k[zm \ m G M^] indexed by m. We make two observations: 
Lemma 3.7. Let mi,mj G M^. With the notation as above 

{{) L{mi) = L{mj) if and only if mj = cr*(mi) for some t G {0, . . . , n — 1}. 
(m) //deg(gcd(m, a(m), . . . , cr"~^(m))) = d—1, then the class of L{m) is a zero divisor 
on T/I{Sn{d)). 

Proof, (i) Suppose that L{mi) = L{mj). So 

So Zrrij = Zat{mi) for some integer t, i.e., ruj = a^{m.i). Conversely, if ruj = cr*(mj), then for 
all / G {0, . . . , n — 1}, cr'(mj) = (T'((T*(mj)) = cr'+*(mj) = a^{mi) were p = l+t (mod n — 1) 
since a has order n — 1. The conclusion now follows. 

(ii) We begin with the observation that if then 

_ ^'*i+V"!+2 . . . ai . . a„ 

Thus, if a = min{ai, . . . , a„_i}, then gcd(mi, cr(mi), . . . , a^-'^^nii)) = x'J ■ ■ ■ x'^_^x'^". This 
monomial will only have degree ci — 1 if n — 2 of {ai, . . . , a„_i} equal a, and the remaining 
exponent equals a + 1. Say that Oj = a + 1. 

Consider the monomial rh = = xl ■ ■ ■ x'^_ix'^'^^^ . Let Zm be the corresponding 
monomial in T. The linear forms L{m) and Zfh are not in I{Sn{d)). However, 

ZrhL{m) = Zfh{z.m + Z„{^rn) H h Z^n^H^rn)) G I{Sn{d)) 

because m and cr*(m) are adjacent in Sn{d) for each z = 0, . . . , n — 1. This implies {ii). □ 

Our second strategy for constructing C is to form the set 

C = {L{m) I m G Ma} \ {L{m) \ deg(gcd(m, a(m), . . . , (T"~-^(m))) = d - 1}. 

By Lemma [3771 (z) we can pick all the L(m) to be distinct. Moreover, we want to eliminate 
those linear forms that are zero divisors since they will never be part of partial hsop. 

In contrast to our previous method, we found that the computation of dimT/ J with J = 
I{Sn{d)) + (L I L G £) was significantly easier than the computation of dimT/ I{Sn{d)). 
It also required less memory. However, this method has the disadvantage of introducing 
too many extraneous linear forms, resulting in large upper bounds. 
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3.3. Neighbours. Our third choice for constructing £ makes use of the neighbours of a 
vertex. Recall that for any m G M^, the neighbourhood of m, is the set N{m) = {m' G 
Md I {m,m'} is an edge of Sn{d)}. To each m G we associate the linear form 

L{m) = Zm+ ^ Zra'. 

m'&N{m) 

Computer experimentation has then suggested the following construction of C Suppose 
that some oracle has given some maximal independent set W of vertices in the graph Sn{d) 
(in the next section, we describe an algorithm which does this). Our desired set is 

C = {L{m) : G W]. 

When we applied Algorithm 13.51 with this set of linear forms in Step 2, we found that 
for known values of n and d (i.e., small n and d), our bounds were close to the true 
value of an{d). In other words, C was almost a partial hsop. At the same time, the 
computation of dimT/J with J = I{Sn{d)) + {L \ L E C) appeared easier than the 
computation of dimT/ I{Sn{d)). Unfortunately, as n and d increase, the computation of 
dimT/J bececomes difficult; see Section 5 for more details. As a final comment, it would 
be of interest to find a good theoretical reason why one should use a maximal independent 
set. 

4. Greedy Algorithms for Computing Bounds 

In this section, we use the symmetry of the graph Sn{d) to either compute or place 
bounds on a„((i) and Pn{d)- The symmetry is used to reduce the number of cases we 
must check in our algorithms. We first present a greedy lower bound algorithm for an{d). 
Following that, we will explain some additional ways to use the symmetry of Sn{d) and 
present two greedy upper bound algorithms for pn{d + 1). 

4.1. The Spreading Number. Be Lemma [2.41 a„ (d) is the cardinality of the maximum 
independent set on Sn{d). Finding a random maximal independent set is trivial with a 
greedy algorithm. Depending upon our initial choice of vertex and the graph itself, the 
difference between the cardinality of the resulting set and an{d) will vary greatly. The 
following algorithm is well known; we include it for completeness. 

Algorithm 4.1. Compute a lower bound for a„((i) 

Input: n,d — The number of variables and degree of monomials, respectively. 

W — An initial independent set. (Optional) 

Output: The cardinality of a maximal independent set on Sn{d) 

Step 1 Let L denote the set of all vertices not in W and not adjacent to any m G W. 

Step 2 Select the vertex v E L with the smallest degree, and add it to W. 

Step 3 Remove v and N{v) from L. If L = 0, return \W\. Otherwise, return to Step[2j 

We describe an alternative approach to finding a lower bound for a„((i). The moti- 
vation comes from an observation about the case n = 3. Interior vertices of 5*3 (c?) are 
of degree 6, and if we regard the vertices of a maximum independent set as the centres 
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of hexagons, then a maximum independent set on S^ld) yields a hexagonal tiling of the 
plane. Constructing this sort of tiling on 5*3 (ci) therefore yields a lower bound for 03 (ci). 
In the next algorithm, we pick a vertex "in the middle" of Sn{d). We then work out 
towards the exterior of Sn{d)] the vertices in our independent set can be thought of as 
the centre of ra- dimensional tile. 

Algorithm 4.2. Compute a lower bound for a„((i) 

Input: n,d — The number of variables and degree of monomials, respectively. 
Output: The cardinality of a maximal independent set on Sn{d) 

Step 1 Let = and E = M^. 

Step 2 Let a = [^J , and set f 1 = X1X2 ■ ■ ■ x^.^x^ with e = d — {n — l)a. Designate Vi the 
origin. Add vi to W. 

Step 3 Let E := E \ (A(t>i) U {fi}) be the set of all vertices eligible for inclusion in 
W. Choose a vertex V2 one degree removed from the origin, i.e, adjacent to some 
element of N{vi). Then V2 is the new origin. 

Step 4 Repeat the above step. At each step i, choose the new origin to be the vertex with 
the most neighbours already discarded from E. This preserves as many vertices as 
possible for consideration. If there are no eligible vertices one degree removed from 
the current origin, backtrack to the previous origin, and continue until E = ^. 

Step 5 Return \W\. 

4.2. The Covering Number. By Lemma [2.41 p„ id 4- 1) is the minimum upward clique 
cover of S'„(c/+1). As with finding maximal independent sets, we can use greedy algorithms 
to produce minimal clique covers of Sn{d+1). These provide an upper bound for p„((i+l). 

By definition, every upward clique is uniquely identified with a monomial from M^. For 
vertices of Sn{d + 1) that consist of more than one indeterminate, many factorizations 
are possible; e.g., Xix\ can be written as {xi)x\ or {xiX2)x2- For monomials of the form 
xf"*"^, however, there is only one such factorization: xfxi. Thus, x^"*"^ belongs only to the 
upward clique identified with xf . These unique upward cliques containing each xf"^^ must 
therefore by in any clique cover of S'„((i + 1), so we can use them as the initial set when 
constructing an upward clique cover. 

Algorithm 4.3. Compute an upper bound for Pn{d). 

Input: n,d — The number of variables and degree of monomials, respectively. 
Output: A minimal upward clique cover of Sn{d) 

Step 1 Initialize our cover C with the set of upward cliques that contain xf. 
Step 2 Select an uncovered vertex, v, at random. If there are no vertices left uncovered, 
go to Step m 

Step 3 Iterate over the upward cliques that contain v. Select the upward clique with the 
fewest number of vertices already covered and add it to C. Return to Step [2J 

Step 4 For each v G M^, compute its frequency, i.e., the number of upward cliques that 
contain it, in C. 

Step 5 Iterate over the elements of the C. If an upward clique does not contain a vertex 
of frequency 1 — all its vertices are represented by other cliques as well — then it 
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is not essential to the cover, so discard it. Repeat this step until we complete an 
iteration without discarding any cliques. 
Step 6 Return C minimal cover. 

Aside from our choice of initial members of the cover, this method does not take into 
account any specific information about Sn{d+1). Curtis [2J proved a formula for p^^d + l) 
using combinatorial methods, and |4] has a proof of the same result using linear program- 
ming. Hulett and Will also provided lower bounds for p^^d + 1), when d > 25, and 
upper bounds for p4^{d + 1), when d > 4. They observed that, owing to the structure of 
Sn{d), vertices on the interior of the graph can be covered exactly once, and the "excess" 
members of the cover come from vertices along the outside of the graph. 

We adapted our greedy algorithm to this idea, removing the random selection of vertices 
to cover and instead working outward, covering the interior vertices first and the vertices 
closest to the outside last. We also use Sym(?T,) to create orbits of the vertices of Sn{d). 

Definition 4.4. For any m G M^, the orbit of m is the set 

{cr(m) I cr G Sym(n)}. 

Since elements of Sym(n) do not alter the exponents of a monomial, only the order of the 
exponents relative to the indeterminates, the orbit of m is also the set of all permutations 
of the exponents of m. By definition, the exponents of any m G Md always sum to d, 
and therefore the orbits of Sn{d) are in an one-to-one correspondence with the integer 
partitions of d of length at most n. We can write orbits as vectors in N", and in this 
form it is easy to determine whether an orbit is an independent set, a clique, or neither 
by examining the entries in the vector. Iterating over the list of orbits of Sn{d + 1) in 
reverse lexicographical order will help us in computing an upper bound on p„((i + 1). 

Algorithm 4.5. Compute an upper bound for pn{d). 

Input: n,d — The number of variables and degree of monomials, respectively. 
Output: A minimal upward clique cover of Sn{d) 

Step 1 Initialize our cover C with the set of upward cliques that contain xf. 

Step 2 Obtain a list, L, of the orbits of Sn{d), where each orbit is represented as a vector 

in N". Sort the list in reverse lexicographical order. 
Step 3 Iterate over L. For each orbit O & L, iterate over the vertices v E O. If f is 

covered, continue. If not, iterate over the upward cliques containing v. Select the 

upward clique that contains the fewest number of vertices already covered, and 

add it to C. 

Step 4 Minimize the cover by computing the frequency of each vertex and eliminating 

unnecessary upward cliques, as described in Algorithm 14.31 
Step 5 Return C clS 8b minimal cover. 

4.3. From Bounds to Exact Values. The greedy algorithms described for both a;„((i) 
and Pn{d+ 1) can be extended, using a brute force approach, to compute the exact value 
of these numbers. For a set W of size u < an{d), we can iterate over the sets of size u + 1 
to check if a larger independent set exists. If we do not find such a set, then obviously 
C(n{d) = u; otherwise, we can maximize the new independent set and repeat the process. 
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Similarly, for a cover C of size t > p„((i + 1), we can iterate over families of cliques of size 
t — 1 to check for a smaller cover. With each iteration, the number of combinations we are 
required to check grows quite quickly. We can reduce this number using the symmetry of 
the graph: each set T C is associated, via Sym(r;,), with n\ (not necessarily distinct) 
other subsets of M^. Hence, when we eliminate T, we also eliminate up to n\ other 
potential sets that we must check. In practise, this improvement is not significant enough 
to make this method feasible. For example, the best lower bound we have for as (5) is 26. 
Since IM5I = 126, we would have to check {^^^^ ^ 2.33 x 10^^ sets to ascertain whether 
o:n{d) = 26. This is an improvement of several orders of magnitude upon the purely 
brute-force approach, where we would need to check 8.5 x 10'^^ sets. However, it is still 
computationally infeasible when run as a serial process, even on an HPC cluster provided 
by SHARCNET. 



5. Results and Comparison to Known Bounds 

In this section we will summarize the bounds for a;„((i) and pn{d+l) reported in previous 
papers and compare them to bounds we found with our algorithms from Sections [3] and HI 
We conclude the section with two tables that present the known values of the spreading 
and covering numbers, as well as the best known bounds for certain n and d. Some of the 
computations were performed in Macaulay2 1.3.1 with 4 GB of memory allocated to 1 
CPU and 1 node on SHARCNET's Saw cluster!! Other computations ran on the Kraken 
clusteiH and used up to 16 GB of memory in Macaulay2 1.4. The run times are taken 
from Macaulay2^s time function. 

5.1. Existing Bounds. As well as providing the first definition of the spreading and 
covering numbers, Geramita, Gregory, and Roberts proved the following inequalities: 

Theorem 5.1 ([S^, Theorem 5.2 and Proposition 5.9]). For all n > 2,d > 2, we have 

(1) ^ < Md) < < ^ + ^^v^.^id), and 

(2) ^-^<Pn(^)- 

Moreover, they presented formulas for a^i^d) |3j, Theorem 5.4], Q;n(3) and Pn(3) [31 
Theorem 5.8], and a^^i^d) for d odd [3, Proposition 5.6]. They include a formula, without 
proof, for the case of even d as well. We will refer to these formulas in Section [6] to 
prove the correspondence of a^^d) to an interesting integer sequence. Curtis [2J proved a 
formula for p^{d + 1), and in [3], Hulett and Will presented an alternative proof of this 
formula, as well as bounds on pi{d+ 1): 

Theorem 5.2 ([4^ Theorems 4.1 and 4.2]). For all d > 5, we have 

(1) ifd IS odd, pi{d) < {d^ + 15d'^ - 6ld + 261)/24, or 

(2) ifd is even, pi{d) < [d^ + 15d^ - 3Ad + 240)/24. 



https : //www. sharcnet . ca/iny/systems/show/41 



https : //www . sharcnet . ca/iny/systeins/show/69 
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5.2. Comparison of Bounds for an{d). We compared the lower bounds of [3] to our 
computations witli Algoritlim 14. 2[ as well as our random greedy algorithm. Like an{d) 
itself, the spreading number's existing upper bounds are similarly difficult to compute, 
so we focused on how closely Algorithm 13.51 approximates known values of In 
Tables [T] and m GGR refers to the lower bound for an{d) from in Theorem I5.1IIT] ) , and 
12.71 refers to the recursive lower bound from Theorem 12.71 while 14.11 and 14.21 refer to 
the bounds computed by those respective algorithms. Recall that Algorithm 14.11 has a 
random component, and so the reported result is the maximum taken over 100 iterations. 



d 


asid) 


GGR 


12.71 


14.11 


14.21 




d 


0:4 ((i) 


GGR 


12.71 


14.11 


14.21 


3 


4 


4 


3 


4 


4 




3 


5 


5 


5 


5 


5 


4 


6 


5 


6 


6 


5 




4 


11 


9 


10 


11 


9 


5 


7 


7 


7 


7 


7 




5 


14 


14 


12 


14 


13 


6 


10 


10 


10 


10 


10 




6 


24 


21 


21 


24 


19 


7 


12 


12 


11 


12 


12 




7 


30 


30 


26 


27 


27 


8 


15 


15 


15 


15 


15 




8 


45 


42 


39 


40 


35 


9 


19 


19 


17 


19 


19 




9 


55 


55 


49 


48 


48 


10 


22 


22 


21 


20 


22 




10 


76 


72 


67 


64 


60 



Table 1. Comparison of Table 2. Comparison of 

lower bounds for a^^d). lower bounds for a^i^d). 



The GGR lower bound remains the best lower bound. For n = 3, the two algorithms 
perform quite well compared to this bound, while the bound of Theorem 12.71 is noticeably 
worse. For n = 4, the algorithms remain comparable with each other but are still worse 
than the GGR bound. Algorithm 14. 2[ which begins at the centre of the graph and works 
outward, produces worse lower bounds than the greedy Algorithm 14.11 Both algorithms 
are similar in terms of memory consumption and time complexity. As an example of the 
time needed. Algorithm 14.11 takes 0.40 seconds to compute a bound for a4(3) but 40.99 
seconds to compute a bound for a4(10). 

Overall, we have not made much progress improving lower bounds for Our 
algorithms are useful if one needs to generate a maximal independent set on Sn{d), but if 
one needs only a bound, then the GGR formula remains superior. 

Comparing our upper bounds from Section |3] with the bounds from [3] is more difficult. 
Geramita, Gregory, and Roberts prove several upper bounds that, like an{d), are linked 
to properties of Sn{d) that make them difficult to compute. Our algorithms avoid this 
by using commutative algebra techniques; unfortunately, these methods introduce their 
own performance issues. However, our algorithms can compute an upper bound that, in 
tandem with a lower bound, provide a specific range into which each of must fall. 

We tested two versions of Algorithm 13.51 differing in our choice of linear form. The first 
version uses the symmetry of Sn{d), as described in Section |3^ The second version uses 
neighbours of a vertex, as described in Section 13.31 We encountered practical limitations 
with both versions. The first version does not produce very tight bounds; for example, 
for 04(8) = 40, the algorithm gives 79 as an upper bound. The second version produces 
very good bounds, at least for the few values we could compute. This version consumes 
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a great deal of memory. For n = 3,3 < d < 9 this version gives bounds that match the 
known value of asld). We ran out of memory computing a bound for 0:3(10), and the 
same is true for n = 4, (i > 4, and n = 5,Q,d > 3. 

5.3. Comparison of Bounds for pn{d + 1). We will examine how the upper bounds 
found in [3j and [4J compare to our computations with Algorithms 14.31 and 14.51 In Table 
El GGR refers to the upper bound for pn{d+ 1) from ([T]) in Theorem 15 . 1 1 and HW refers 
to the bounds from Theorem 15. 2[ while 14.31 and 14.51 refer to the bounds computed by 
those respective algorithms. 



d 


GGR 


HW 


14.31 


14.51 


5 


30 


19 


22 


19 


6 


42 


33 


33 


29 


7 


57 


38 


44 


40 


8 


75 


60 


62 


55 


9 


97 


69 


81 


74 


10 


121 


100 


113 


96 


11 


150 


114 




122 


12 


182 


155 




147 


13 


219 


175 




185 


14 


260 


227 




223 


15 


306 


254 




275 



Table 3. Comparison of upper bounds for pi{d). 

All three revised bounds are better than the original GGR bound. Our original greedy 
algorithm, 14. 3[ fails to match HW. However, by taking advantage of the symmetry and 
structure of Sn{d), our revised algorithm. Algorithm 14. 5[ is quite close to HW. In fact, it 
seems that for even d our bounds are equal or better, with the reverse is true for odd d. 
This pattern holds for at least d < 24, with the exception of ci = 22. We are not certain 
why this is the case. Finally, recall that the HW bound holds only for n = 4, whereas our 
algorithm works for all n > 2. When tested against GGR for small values of d for n = 5, 6, 
Algorithm 14.51 consistently performs better (refer to Table [5] for example output). 

Neither of our algorithms consume much memory, but as one might expect, as d in- 
creases the computational time increases significantly. As a result, we did not compute 
bounds for d > 10 with Algorithm 14.31 the computation for d = 3 took only 2.37 seconds, 
but for d = 10, the algorithm took 1709.44 seconds to terminate. Algorithm 14.51 once 
again performs better, taking 0.05 seconds for d = 3 and 21.28 seconds for d = 10. With 
each increase in degree, the amount of time to compute a bound roughly doubles, with 
d = 24 taking 83051.40 seconds. This complexity increases as n increases; for example, 
we computed a bound for pe{3) in 2.77 seconds, but a bound for P6(4) takes 7.81 seconds. 

Our computations indicate that our algorithms yield tighter upper bounds for p„((i + 1) 
than those of [3] and, n = 4 and d even, [4]. In some cases, it might be preferable to 
sacrifice some accuracy in favour of speed, particularly for values of n and d at which our 
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algorithms might become infeasible. Even so, we hope this provides a useful example of 
how one can use the structure and symmetry of Sn{d) along with a greedy algorithm to 
improve bounds on an{d) and pn{d + 1). 

5.4. Tables of Results. These tables are based on those from [Ij. The a4^{d) row was 
filled according to the formulas in [3] . In some entries we have used a bracketed notation, 
[a,b], to indicate a range given by known lower and upper bounds a and b, respectively. 
For values we could not compute exactly, we have attempted to bound them within a 
neighbourhood [a, b] where |6 — a| < 100. If no value is given, then we could not find such 
a neighbourhood. 



d 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


ai{d) 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


a2{d) 


1 


2 


2 


3 


3 


4 


4 


5 


5 


6 


aaid) 


1 


3 


4 


6 


7 


10 


12 


15 


19 


22 


a4{d) 


1 


4 


5 


11 


14 


24 


30 


45 


55 


76 


a5{d) 


1 


5 


7 


16 


[26,33] 


[42,61] 


[66,94] 


[99,142] 


[143,209] 


[201,285] 


ae{d) 


1 


6 


10 


[24,30] 


[42,63] 


[77,119] 


[132,201] 


[215,] 


[334,] 


[501,] 


aj{d) 


1 


7 


14 
















a%{d) 


1 


8 


[16,20] 
















ag{d) 


1 


9 


















aio{d) 


1 


10 



















Table 4. Computed values and bounds of the spreading number 



d 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Pl{d + 1) 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


P2(d+1) 


2 


2 


3 


3 


4 


4 


5 


5 


6 


6 


P3(d + 1) 


3 


4 


6 


9 


12 


15 


18 


23 


27 


32 


P4(d + 1) 


4 


6 


12 


18 


27 


[32,38] 


[45,55] 


[58,69] 


[77,96] 


[96,114] 


P5(d+1) 


5 


9 


20 


33 


[42,61] 


[66,94] 


[99,142] 


[143,209] 


[200,285] 


[273,] 


Peid + l) 


6 


12 


30 


[42,63] 


[77,119] 


[132,201] 


[214,] 


[333,] 


[500,] 


[728,] 


Prid+l) 


7 


16 


[30,46] 


[66,102] 


[132,210] 












psid+l) 


8 


20 


















P9(d + 1) 


9 


25 


















Pwid+1) 


10 


30 



















Table 5. Computed values and bounds of the covering number 



6. A053307 

We conclude by showing a relationship between a^^d) and a known integer sequence. 
For 0^4 (d), [3] gives an explicit formula; in particular 

fA^ — / ^^^^ ^ even 

Md) - ^ 3^ ^ 

Here, v^^d) is the number of vertices of S^i^d). 
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Theorem 6.1. For d > a4^{d) is the number of non-negative integer 2x2 matrices 
with sum of entries equal to d, under row and column permutations. 

Proof. Let a{d) be the number of non-negative integer 2x2 matrices with sum of entries 
equal to d, under row and column permutations. This sequence is in the OEIS as A053307 
[7]. In fact, it is an interleaved sequence where 

a(2d + 1) = A000330(d + 1) = ('^ + ^K'' + 2)(2d + 3) 

6 

and 

, „ , ^ (d+lY + 2(d+l) 

a{2d) = A006527(rf + 1) = ^ ^ ^ \ 

3 

Recalling that V4{d) = {^~^^), it follows that 

a,i2d+l) = ^'^^'^'^^^ - ^''^'^ - (2rf + 4)(2rf + 3)(2rf + 2) 



4 4 3!(4) 

(rf + 2)(2d + 3)(rf+ 1) 



6 

Similarly, 



a{2d+ 1). 



"^(2rf) - + ^ + ^ 

_ 8d^ + 24^2 + 22rf + 6 6d + 6 _8d^ + 2Ad^ + AOd + 24 
~ 24 ^ 8 ~ 24 

_ d^ + 3d^ + 5d + 3 _ {d+lf + 2{d + l) _ 



□ 



Even though the sequence A053307 and 0:4 (d) are related, it is not immediately apparent 
why this sequence and the spreading numbers of 6*4(6?) are linked. The correspondence 
may be a result of the two interleaved sequences that make up A053307: the square 
pyramidal numbers, for odd d, and the sum of two tetrahedral numbers, for even d. Thus, 
we conclude with the following question. 

Question 6.2. What is the relationship between the spreading number in four variables 
and the sequence A053307? Why are they equivalent? 

Hopefully explaining the relationship between 04 (rf) and A053007 could open up new 
techniques for computing the spreading and covering numbers. 
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